Swimming active droplet: A theoretical analysis 
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PACS 47 . 55 . D- - Drops and bubbles 
PACS 47 . 55 . pf - Marangoni convection 

Abstract - Recently, an active microswimmer was constructed where a micron-sized droplet of 
bromine water was placed into a surfactant-laden oil phase. Due to a bromination reaction of the 
surfactant at the interface, the surface tension locally increases and becomes non-uniform. This 
drives a Marangoni flow which propels the squirming droplet forward. We develop a diffusion- 
advection-reaction equation for the order parameter of the surfactant mixture at the droplet 
interface using a mixing free energy. Numerical solutions reveal a stable swimming regime above 
a critical Marangoni number M but also stopping and oscillating states when M is increased 
further. The swimming droplet is identified as a pusher whereas in the oscillating state it oscillates 
between being a puller and a pusher. 
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Introduction. — A rigorous understanding of swim- 
ming on the micron scale is crucial for developing microflu- 
idic devices such as a lab-on-a-chip [T]. This understand- 
ing comes from watching nature, i.e., by studying the loco- 
motion of living organisms such as bacteria or algae [2] but 
also from designing artificial microswimmers, used for ex- 
ample as medical microrobots [3] . Both, real live cells and 
man-made microswimmers, have thoroughly been used to 
study interaction between swimmers [4], interaction with 
walls [SHS], or swarming [9]. One possible design of an 
artificial swimmer is an active droplet. Here, we think 
of a droplet with a surface where a chemical reaction 
occurs. Alternatively, droplets or bubbles can be made 
active by having an internal heat source [10| . Droplets 
are particularly interesting systems to study since they 
are used extensively in microfluidic devices as microreac- 
tors in which chemical or biological reactions take place 
[TTJ[T2]. In the following we give an example of an active 
droplet and investigate in detail its propulsion mechanism. 
Self-propelled active droplets have been studied in vari- 
ous experiments, including droplets on interfaces [TS l fM j . 
droplets coupled to a chemical wave [15] , and droplets in a 
bulk fluid [T6HI9] . Theoretical treatments include a model 
of droplet motion in a chemically reacting fluid [5D] , stud- 
ies of the stability of a resting droplet [211 - |25| , and simu- 
lations of contractile droplets and of droplets driven 
by nonlinear chemical kinetics |27j . 

The swimming active droplet we consider in the follow- 
ing is a solution of water and bromine which is placed in a 



surfactant-rich oil phase [TH] . The resulting water droplet 
has a typical radius of 80^m. In order to lower the sur- 
face tension and thus the total energy of the system, the 
surfactants in the oil phase form a dense monolayer at the 
droplet interface, giving the droplet the structure of an 
inverse micelle. The observed directed swimming motion 
of the droplet with a typical swimming speed of 15/im/s 
can be understood as follows [TS|. 

The bromine within the droplet chemically reacts with 
the surfactants in the interface which results in a weaker 
surfactant. Hence, the 'bromination' reaction locally leads 
to a higher surface tension in the interface. As a conse- 
quence local gradients in surface tension will lead to a fluid 
flow at the interface and in the adjacent fluid inside and 
outside of the droplet in the direction of increasing sur- 
face tension. This effect is called Marangoni effect. The 
fluid flow then leads in turn to advection of surfactants at 
the interface. As a result gradients in surface tension are 
enhanced. Thus, the resting state becomes unstable and 
the droplet starts to move. Additionally, brominated sur- 
factants arc constantly replaced by non-brominated sur- 
factants from the oil phase by means of desorption and 
adsorption. The droplet stops to swim when cither the 
bromine or the non-brominated surfactants in the oil phase 
are exhausted. This was also observed in the experiments 

m- 

The active droplet is an interesting realization of the 
'squirmer' [281129] which has been introduced to model the 
locomotion of microorganisms. Often they propel them- 



p-1 



M. Schmitt et al. 



selves by a carpet of beating short filaments called cilia on 
their surfaces. Instead of modeling each cilium separately, 
one prescribes the fluid flow at the surface initiated by 
the beating cilia which then drags the squirmer through 
the fluid. Here, for the active droplet the surface flow is 
generated by the Marangoni effect. 

The swimming active droplet crucially depends on the 
dynamics of the mixture of non-brominated and bromi- 
nated surfactants at the interface. In this article we 
model it by means of a diffusion-advection-reaction equa- 
tion based on a free energy functional for the surfactant 
mixture. Numerical solutions then show that in a certain 
parameter range the resting state of the droplet becomes 
unstable and the droplet starts to move. The solutions 
reach a stationary state corresponding to a swimming mo- 
tion and conflrm that the droplet is a pusher [5] , as found 
in the experiments [18) . In addition, we identify further 
patterns of motion. We find that the droplet stops after 
an initial motion or that it oscillates back and forth. 

Model. — In order to model the droplet propulsion we 
set up a dynamic equation for the surfactant mixture at 
the droplet interface that includes all processes mentioned 
before. We assume that the surfactant completely cov- 
ers the droplet interface without any intervening solvent. 
We also assume that the area of both types of surfac- 
tant molecules (brominated and non-brominated) is the 
same. Denoting the brominated surfactant density by Ci 
and the non-brominated density by C2, we can therefore 
set ci + C2 = 1. We then take the concentration difference 
between brominated and non-brominated surfactants as 
an order parameter (j) = ci — C2. In other words 0=1 
corresponds to fully brominated and </> = — 1 to fully non- 
brominated and ci = (l4-0)/2andc2 = (1 — (^)/2. Finally, 
we choose a constant droplet radius R. 

Dijfusion-Advection-Reaction equation. The dynam- 
ics of (f) is governed by a diffusion-advection-reaction equa- 
tion: 

4' ^ -V ■ (in + U) - Tj^.\(f' ~ 4'eq) , (1) 

with diffusive current j^) and advective current due 
to the Marangoni effect. The third term on the right- 
hand side of Eq. ^ is the reaction term and describes the 
bromination reaction as well as desorption of brominated 
and adsorption of non-brominated surfactants to and from 
the outer fluid, tr is the timescale on which these pro- 
cesses happen and t^eg sets the equilibrium coverage of 0. 
In other words, ad- and desorption dominates for 0eg < 0, 
while bromination dominates for (p^q > 0. Imagine for in- 
stance the case (jjeq ~ 1, i.e., a droplet with bromination 
but without ad- and desorption of surfactants. The reac- 
tion term would then always be positive, therefore driving 
the droplet to a completely brominated state 0=1. 

The general mechanism of Eq. ^ is as follows. The dif- 
fusive current always points 'downhill', oc — V0. How- 
ever, we will show below that the opposite is true for 
since approximately cx V0. Thus, apart from the reac- 
tion term, jjy and are competing and as soon as 



dominates over j^i, experiences 'uphill' diffusion, i.e. 
phase separation. As a result the resting state will be- 
come instable and the droplet will start to move. We will 
now present a careful derivation of ju and from a free 
energy approach. This shows that the diffusive and advec- 
tive currents in Eq. (jlj are in general non-linear functions 
of 0. 

Diffusive current. The basis for the following is a free 
energy density / for the droplet interface, which we write 
down as a function of concentrations ci and C2. In for- 
mulating the free energy density /, we follow the Flory- 
Huggins approach |30] . Accordingly, / is composed of the 
mixing entropy plus terms mimicking lateral attractive in- 
teraction between surfactants: 



/ 



A 



[ci In ci + C2 In C2 — bic\ - 62C2 - 612C1C2 



(2) 

where A denotes the area of a surfactant in the inter- 
face and bi (62) is a dimensionless parameter characteriz- 
ing the interaction between brominated (non-brominated) 
surfactants and 612 the interaction between different kind 
of surfactants. In the following we assume for simplicity 
bi2 = {bi + &2)/2. In terms of the order parameter we 
obtain: 
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i±^lni±^ + i^lni^ 



(3) 



-i(&i + fo2)- 1(61-62) -^(61 + 62) . 
The total free energy is then given by the functional 



f{cj>)dA 



(4) 



For a conserved order parameter field the diffusive cur- 
rent is proportional to the gradient of the variation in F 
with respect to [31] : 



-AV^ = -A/"(0)V0 
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(5) 

(6) 



with positive mobility A. Substituting j/^ into Eq. ([T]) 
yields a Cahn-Hilliard type equation [32]. Note that the 
diffusion constant in Eq. (|6]) decreases with increasing in- 
teraction energy. In fact, the condition ju oc — V0 is only 
fulfilled for a convex free energy with /"(0) > 0, i.e. if 
61 -I- 62 < 4. In addition, the diffusion coefficient in j^i 
is smallest for = 0. It increases with |0| and diverges 
at 101 = 1. An alternative approach of deriving diffusion 
currents in mixtures is presented in |33[l34j. 

Advective current. The advective current for the order 
parameter is given by 



JA 



(7) 
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Swimming active droplet 



where u is the velocity of the surfactants at the interfaced 
Since we are studying the active droplet in an axisym- 
metric geometry, we assume = (f>{0) and u = ug{9)eg, 
where the front of the droplet is at = 0, see inset of 
Fig. [T] (b) . For this geometry there exists a solution of the 
Stokes equation for the fluid flow field inside and outside 
of the droplet as well as the fluid velocity at the interface 
[35j|36]. The solution at the interface is given in terms of 
the surface tension gradient: 



E 

Tl = 2 



n{n — 1) 
2ry 



where z = cos(0). r/ = rii + r^o is the sum of the viscosities 

— 1/2 

inside and outside of the droplet and Cn arc Gcgen- 
bauer polynomials. They are related to Lcgendrc polyno- 
mials by Pn{z) = —j^C~^{^{z). Equation ^ is nothing 
but a representation of the the Marangoni effect. It essen- 
tially states u oc Vcr, i.e., a fluid flow in the direction of 

V(T. 

Thus, in order to calculate ue, we need an expression 
for da/d9, which can be found by deriving an equation 
of state for the surface tension a. The surface tension a 
is the thermodynamic force conjugate to the surface area. 
This gives: 

f)f rif 

(9) 



, Of df 

OCi OC2 



which we identify as the Lcgcndre transform of the free 
energy ([2]) to the chemical potentials fii = Hence, 
a = [bicf + 62C2 + bi2CiC2] , or in terms of cj) and 
again with 612 = (fei + &2)/2: 



l{bi + 62) + 2(61 - b2)(f> + lib, + b2)cf>^ 
o o 



(10) 

In order to obtain the proper behavior of the equation 
of state, i.e. an increasing surface tension with increas- 
ing (j>, we need to assure that cr'(^) > 0. This holds 
if hi > 62, meaning that the interaction energy between 
brominated surfactants has to be higher than between the 
non-brominated ones. Note that in the limit of — the 
equation of state becomes linear in 0. The gradient of a 
is given by 



2A 



(bi - 62) 



1 + - 



7 61 + 62 



de 



8 61 - 62 

(11) 

By substituting this into Eq. one can calculate the 
advective current ([7|) for a given (t>{0). 

Eqs. ^ and PT|) essentially state that u oc V0. There- 
fore, when > 0, the advective current = 4''^ appar- 
ently always points 'uphill', i.e., in the opposite direction 
compared to jo- On the other hand, when < 0, the 



^Let the advective currents of the two types of surfactants be 
= ciui and = C2U2. Under the assumption that the individual 
velocities are identical ui = U2 = u, one obtains 
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Fig. 1: (a) Stationary order parameter profiles after 10^ time 
steps for (paq = 0.5 and several Marangoni numbers M. Gray 
solid line: Initial condition, (b) Corresponding interface veloc- 
ity profiles. Inset: Droplet geometry. 



advective current acts 'downhill', i.e., in the same direc- 
tion as jD. As a consequence, the Marangoni flow will 
only drive the droplet when (f> > 0. This is the case 
when there are more brominated surfactants than non- 
brominated ones. 

Together with ^ and ([7]), Eq. ^ becomes a closed 
equation for (j). Writing gradients in units of and time 
in units of the diffusion time scale td = A{\kBT)~^ 
yields 



= -V • [io + M(t>VL) - k{ 



'■'eq 



(12) 



where the currents j/3 
less and 



M = 



and jyi = M(j)u are now dimension- 
(61 - b2)R 



Xr, 



(13) 



is called Marangoni number. This number compares the 
advective current due to the Marangoni effect, which scales 
as ksTibi — b2){RAri)^^ , to the diffusive current. Ac- 
cordingly, K = T£)Tjj^ is the ratio between diffusion and 
reaction time scale. 

Results. — Wc numerically solve the diffusion- 
advection-reaction equation for with the initial condi- 
tion 0(6*) = (peg + S(l){9), where S(t){d) is a small perturba- 
tion [solid line in Fig. [TJa)]. The boundary conditions at 
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Fig. 3: Phase diagram of the active droplet in {(j>eq, M) param- 
eter space. Examples for the order parameter profiles at the 
positions marked with numbers are given in Fig. [ija) (regime 
1-3) and Fig. (regime 4). 
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Fig. 2: (a) Droplet swimming velocity vd for swimming, stop- 
ping, and oscillating droplets. Parameters are the same as in 
Fig. [U and case 4 belongs to M — 10.5. (b) Depiction of the 
chemical wave of case 4 in a (j){6,t) plot. Same timescale as in 
(a). 



= 0, TT are given by a vanishing current, ]d + M(j)u ~ 0. 
We keep k fixed to a value of 0.1 for all numerical solu- 
tions and comment later on the impact of k on the results. 
Therefore, we are left vi^itli the Marangoni number M and 
(peq as the crucial parameters. To assure a convex free 
energy, we set 6i + 62 = 3. 

Order parameter and velocity profiles. Figure [ija) 
shows examples of the stationary order parameter pro- 
file for (peg — 0.5 and different values of M together with 
the corresponding interface velocity profiles in Fig. [Ijb) . 
Starting with a small Marangoni number of M = 2.5, the 
order parameter relaxes into the homogeneous trivial so- 
lution (j) = (peg of Eq. ([T|), thus the droplet rests. Above 
a critical Marangoni number, the order parameter evolves 
to a stationary inhomogeneous profile, as Fig. [T] shows for 
M = 3. In parallel, the droplet velocity vd depicted in 
Fig.[2lja) saturates on a non-zero value. The droplet swim- 

sin^(6')/2, is determined 



ming speed is given by vd 



Since Cr 



-1/2 



(cos(0)) 



by the first coefficient of the sum in ([8]) and thus vd 
- sin 9u0d9. Note that in our approach 



8 Vi+Vc 



reaches a stationary value without having to introduce a 
'backward' Marangoni stress, as suggested in [18]. Further 
increasing the Marangoni number to M = 4.5, the droplet 
starts to swim but then stops rapidly. The stationary or- 
der parameter profile becomes symmetric around 9 ~ ■k/2 



and swimming is not possible. Finally, the droplet reaches 
an oscillating state for even higher Marangoni numbers 
where it oscillates back and forth as the droplet swim- 
ming speed in Fig. [D^a) demonstrates. In this case the 
order parameter 4>{9,t) resembles a chemical wave that 
travels back and forward between 6* = and 9 = tt. The 
wave is depicted in Fig. [2][b). Note that the frequency of 
the oscillation increases with M. Finally, we remark that 
from comparing Figs. [T] (a) and (b), it is now apparent 
that indeed Eq. ([8]) essentially gives u oc \7(j). 

Phase diagram. Figure [3] shows the phase diagram 
in {(j)eq,M) parameter space with the four regimes of the 
droplet dynamics: resting, swimming, stopping, and os- 
cillating. Since there is no swimming motion possible for 
negative (peg, as discussed before, we only show the phase 
diagram in the range < (f>eq < 0.80 We find similar 
phase diagrams for smaller values of k. For k = 0.01 
the swimming region increases in size and then shrinks 
again for k = 0.001 until for k = swimming solutions 
are no longer possible. The critical Marangoni number at 
the onset of the swimming regime stays, however, roughly 
constant. On the other hand, for k = 1 and 10, i.e., in the 
limit of fast bromination reaction and exchange of surfac- 
tants, only resting, stopping and oscillating solutions but 
no stable swimming solutions were found. 

Reduced phase space. Due to the axisymmetric geom- 
etry we decompose the order parameter (p into Legendre 
modes 



= ^^«(cos(0))(/)„(t) 



(14) 



ri=0 



^Due to the </> dependent diffusion coefficient in Eq. JSJ, numerics 
requires a much finer grid above (f>eq = 0.8. However, in several tests 
for different values of M no swimming solutions were found above 



: 0.8. 
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Fig. 4: Droplet dynamics in the reduced phase space {4>i,4>2)- 
Black dots show the fixed points for different values of M; 
from dot to dot M increases by 0.2. The red and the green 
line show, respectively, trajectories in the swimming (M = 3) 
and stopping (M = 4.5) state. Inset (a): limit cycle in the 
oscillating state {M = 10.5). Inset (b): map for the swimmer 
type in (0i,(/)2) space classified by the stirring parameter j3 = 
—(f>2/\4>i \ (see main text). 



Together with Eqs. (fTTj) and ([8|) one obtains an expression 
for uq as a function of the mode amplitudes (pi deter- 
mines the swimming speed and 0„>i corresponds to the 
higher modes of ug. In the following, we use the initial 
condition = 0) = (peq ■ In order to investigate the four 
regimes of the droplet dynamics, we plot in Fig.|4]the fixed 
points in (^1,^2) space for increasing Marangoni number 
M at (t)eq = 0.5. For the cases M = 3 and M = 4.5 the 
full trajectories are shown. Note that this illustration is a 
projection onto only two modes of infinitely many modes 
that make up the full phase space of cj). Starting with the 
resting state, one has a stable fixed point at </>! = 0, (/)2 = 
for M < 2.7. Via a subcritical bifurcation the droplet en- 
ters the swimming state at the critical Marangoni number 
M = 2.8. Figure m demonstrates that for the chosen initial 
condition (</>!, (/)2) ~ (0;0) both modes 4>i and 4>2 develop 
non-zero amplitudes at the same critical Marangoni num- 
ber. The trajectory in the swimming state does increase 
its size with increasing M, whereas the swimming speed 
decreases until the droplet reaches the stopping state at 
M = 4.3. As already observed in Fig. [T] the second mode 
02, which is symmetric around = 7r/2, clearly domi- 
nates in the stopping state. In the oscillating regime above 
M = 9.5 a stationary solution does not exist. Instead, the 
dynamics follows a stable limit cycle as the inset (a) in 
Fig. m demonstrates for M = 10.5. Finally we remark, 
since the bifurcation is subcritical, the critical Marangoni 
number for the onset of the swimming state depends on 
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Fig. 5: Displacement of oscillating droplet plotted versus time. 
Color of line shows the value of stirring parameter /3. 



the chosen initial condition. For example, starting the nu- 
merical solution at 4>i ~ 4>2 = —0.1 the critical Marangoni 
number is M = 1.7. 

The active droplet as pusher. To describe the ba- 
sic features of a squirming swimmer, it is sufficient to 
only study the first two modes of its surface velocity field 
[2l[l8l[29l|37l[38|. While the the first mode determines 
the swimming velocity, the dimensionless 'stirring' param- 
eter /? = —02/1011 characterizes the swimmer type. When 
/3 is negative, the flow around the droplet is similar to 
the flow around a swimming bacterium such as E. coli. 
Such a swimmer is called a 'pusher' since it pushes fluid 
away from itself at the front and at the back. Accord- 
ingly, a swimmer with /? > is called a 'puller'. The 
algae Chlamydomonas is an example for a puller since it 
swims by pulling liquid towards itself at the front and at 
the back [S]. For (3 — > ±00 the droplet becomes a 'shaker', 
i.e., a droplet that shakes the adjacent fluid but does not 
swim. If /3 = 0, the first mode dominates and propels the 
droplet, as is the case for Volvox algae [9]. The classifica- 
tion of the swimmers according to the 'stirring' parameter 
/3 is illustrated in the inset (b) of Fig. Hydrodynamic 
interactions between swimmers and with bounding walls 
depend on their type ('stirring' parameter j3) and strongly 
influence their (collective) dynamics [7l|39]. For instance, 
adjacent pushers generally tend to align and swim parallel 
to each other, i.e., show a polar velocity correlation [HUBT] . 
In fact this kind of behavior was observed in experiments 
of our active droplets [T8| . It is therefore of great interest 
to determine /3. The swimming droplet with (p^q = 0.5 
is a pusher with /3 ranging from —0.7 for A/ = 2.8 to 
— 1.5 for M = 4.2. Similar values from p = —0.5 up to 
—7 were observed throughout the whole swimming regime 
of the droplet. The stopping droplet is always a shaker 
with (3 = —00. Since the limit cycle of the oscillating 
droplet perambulates all four quadrants of the reduced 
phase space, it oscillates in the swimming direction as well 
as in /3, i.e., between being a pusher and a puller. This is 
demonstrated by the droplet displacement plotted versus 
time in Fig. [5l 
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Conclusions. — We have presented a model for an 
active squirming droplet with a surfactant mixture at its 
interface that drives a Marangoni flow and thereby drags 
the droplet forward. Based on a free energy functional 
for the mixture, we derived a diffusion-advection-reaction 
equation for the mixture order parameter at the droplet 
interface. Relevant parameters are the Marangoni number 
M and the reduced relaxation time with which the 
mixture approaches its equilibrium value (peq by bromi- 
nation or de- and absorption of the surfactants from the 
surrounding. 

As predicted from linear stability analysis in [18] . nu- 
merical solutions of the diffusion-advection-reaction equa- 
tion show that above a critical Marangoni number the 
resting state of the droplet becomes unstable. The order 
parameter develops a non-uniform profile and the droplet 
moves with a constant swimming velocity. This only oc- 
curs when the relaxation time (relative to the diffusion 
time) is sufficiently large. The negative stirring parameter 
j3 identifies the droplet as a pusher in agreement with po- 
lar velocity correlations found in experiments |18| . A full 
parameter study in {(jieq^ M) space also reveals a stopping 
droplet, which is a shaker (/3 — — oo), and an oscillating 
droplet that oscillates between being a puller and a pusher. 
We hope that our work initiates further research on the 
active droplet which constitutes an attractive realization 
of the model swimmer called squirmer. 

^ * * 
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